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Abstract 

We analytically compute the Renyi entropies for the RSOS models, representing a wide class of ex- 
actly solvable models with multicritical conformal points described by unitary minimal models and Z„ 
parafermions. The exact expressions allow for an explicit comparison of the expansions around the critical 
points with the predictions coming from field theory. In this way it is possible to point out the nature of the 
so-called "unusual corrections", clarifying the link with the operator content, the role of the symmetries and 
the boundary conditions. By choosing different boundary conditions, we can single out the ground states 
as well as certain combinations of high energy states. We find that the entanglement spectrum is given by 
operators that are not present in the bulk Hamiltonian, although they belong to the same representation of a 
Virasoro Algebra. In the parafermionic case we observe unexpected logarithmic corrections. 
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I. INTRODUCTION 



Entanglement is the unique feature distinguishing a quantum system from a classical one''". 
While we still lack a fundamental, general definition of what entanglement is, we can characterize 
it well when we consider the mutual entanglement of two complementary components of a system 
in a definite state, the so-called bipartite-entanglement . A popular way to quantify it is given by 
the entanglement entropies. In the recent years, it has become increasingly important to be able 
to compute them, either numerically or analytically . Typically, one considers the ground state |0) 
of a quantum Hamiltonian. Once the system has been divided into two parts A, B, it is possible to 
introduce the reduced density matrix, tracing out one of the two subsystems 

PA = TrB|0)(0| (1) 

and then the Renyi entropies are defined as 

5„ = -^lnTrp^, (2) 

I — a 

We notice that, thanks to the free parameter a, the knowledge of the Renyi entropies is equivalent 
to the knowledge of the full spectrum of the reduced density matrix , whose logarithm is known 
as the entanglement spectrum . A particularly important point is the a — > 1 limit, known as Von 
Neumann entropy 

S = lim Sa = —TipA In pa , (3) 

a— 5>1 

which provides a good quantification for the entanglement in terms of a single number. 

Gapped d + 1-dimensional systems obey the so-called area-law : at the leading order in the 
thermodynamic limit of large subsystem sizes, the entanglement entropy is proportional to the 
area of the boundary separating A and B. In d = 1, such law predicts a saturation to a constant 
of the entanglement entropy when A is composed of large intervals. For d > 1 the area law 
remains true for most gapless systems, with possible logarithmic corrections\ These logarithmic 
contributions are a signature of d = 1-dimensional physics. In fact, exploiting the conformal 
invariance of gapless 1 + 1 models, it is known that the entropy grows logarithmically with the 
length of the A interval i, with a proportionality given by the conformal anomaly '. In'", the sub- 
leading contributions were analyzed and the emergence of unusual corrections was linked to the 
effect of relevant (and irrelevant) operators of the critical theory: 

Sa = ^ (—] \ni + c'„ + bj-'^/" + . . . , (4) 
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where c is the central charge of the CFT, c'^ and ba are non-universal constants and x = A + A is 
the dimension of the operator "responsible" for the correction. Marginal operators act differently 
(renormalizing the central charge) and give rise to logarithmic corrections of the form (log£)~^". 

Similar results have been observed"^' ' close to the critical points, where now the correlation 
length ^, being large but finite, becomes the relevant length scale (that is i ^ and have led to 
the conjectured form 

12 \ a J ao 

where ao is a short distance cutoff, Aa and Ba are again non-universal constants and h can be 
interpreted as the dimension of a relevant operator^ In general, one should not expect the entropy 
to be the same scaling function in i and in ^ : although the coefficients of the leading term have to 
coincide, the operators acting in the two cases can be different and x and h can be different. The 
crossover function between the two scaling regime might be accessible using techniques similar 
to'^'". 

Thus, while the entanglement entropy of a gapped d = 1 system could seem not very interest- 
ing, since it saturates to a constant in the thermodynamic limit, its study close to a critical point 
could shed light on the scaling theory governing the lattice models and its universal features. More- 
over, a recently proposed protocol ' would allow the measurement of the Renyi entropies only 
for gapped systems, thus rendering the theoretical computation of the limiting value (5) amenable 
to cold- atom experimental confirmation. 

In this paper we will focus on the analytic computation of the Renyi entropy in this thermody- 
namic limit for the quantum systems obtained from a class of integrable lattice models known as 
Restricted Solid-on-solid (RSOS) . These models and their structures have inspired the discus- 
sion in the last section of ' , on a general relation between the entanglement entropy of quantum 
(integrable) models and Virasoro characters. Inspired by these considerations, we expand and 
detail the calculation sketched in", and extend it to the parafermionic case. The importance of 
the RSOS models is multifold: first of all, they provided the first lattice realization of the unitary 
conformal models as pointed out in-' --. In a different phase, they also realize parafermionic 
models and thus give access to consistent c > 1 CFT's . Moreover, thank to the rich underlying 

' The factors of 2 difference between (4) and (5) is due to the number of boundaries dividing A and B: 2 for an 
interval, 1 in the gapped phase, where both A and B are taken semi-infinite. 
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mathematical structure, they appeared as a fascinating link between integrable lattice models and 
number theory. 

While entanglement is associated to a quantum state, here we will take advantage of the well- 
known link that allows to derive a quantum Hamiltonian from the row-to-row transfer matrix 
of an integrable classical model. In this way, the classical configuration with the lowest free 
energy corresponds to the ground state of the quantum model. Moreover, by a proper choice of 
the boundary conditions, one can select higher energy configurations, which correspond to the 
lowest quantum state within a given sector of the Hilbert space, and is thus a way to investigate 
the entanglement entropy of states other than the ground state. Without showing explicitly the 
quantum Hamiltonian associated to the RSOS transfer matrix, it is worth saying that it naturally 
arises in the context of loop models . More recently, an explicit realization of these Hamiltonians 
has be obtained from a very different perspective as a chain of interacting non-abelian anyons \ 
Another possible approach, whose terminology we decide to adopt here, is to interpret the RSOS 
models as the lattice realizations of an integrable thermal perturbation of a class of rational CFTs. 
This allows to write the action as 

A = AcPT + A j (fx e(x) (6) 

where e(x) is the operator representing the thermal perturbation and A is the coupling constant 
measuring the distance from criticality. 

The RSOS models being ubiquitous and integrable makes the computation of the Renyi entropy 
not only interesting, but also possible analytically by means of the Corner Transfer Matrix (CTM) 
approach. In fact, the reduced density matrix of an half-interval in the thermodynamic limit can 
be shown to be equal (except for the normalization factor) to the CTM operator-"^'-^ 

Pa = 2^ ^pcTM , 2a = TrpcTM > C^) 

and one can therefore compute the Renyi entropy as 

a 1 

Sa = InZi + lnZ„ . (8) 

a — 1 1 — a 

Even though this procedure looks similar to the replica trick exploited in the conformal case, 
here a can be an arbitrary real (or even complex) parameter, thus avoiding all the subtleties of 
the analytic continuation from a = n E N, necessary to compute, for instance, the Von Neumann 
entropy. Therefore, beyond checking that the conjectured form of Eq.(5) applies for the RSOS, our 
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results provide the umpteenth check to the Cardy-Calabrese formula-*^, both for minimal models, 
where the replica trick introduces operators not present in the original Kac table, and also in a 
systems with central charges greater than unity. For minimal models, we identify the leading 
unusual correction of (5) as coming from the second most relevant operator in the model, that 
is A3 3, since the most relevant one, A2,2, is odd under the Z2 symmetry of the ground state. 
However, we find that, by varying the boundary conditions, different sectors can be traced out and 
the leading correction to the entropy can come from other operators as well. Let us stress that this 
result is different from the correction one could naively expect taking the anomalous dimension 
of the perturbation in (6). For parafermionic models, we find the leading correction to come from 
the first thermal field, which is the most relevant only among the 'Lr-2 neutral fields. In addition, 
some boundary conditions at infinity turn on logarithmic corrections, different from those in (5). 
It would be tempting to interpret these terms as due to a marginal operator in accordance with'"''^^; 
however these corrections are present even when the theory does not seem to support a marginal 
field (which is normally related to the existence of a free boson and present only for certain given 
values of r). Thus, the origin of these terms still needs a full explanation and is probably rooted in 
a choice of boundary conditions which has no conformal counterpart in the continuum limit. 

The paper is organized as follows: in section II we introduce the RSOS models, their phase 
diagrams and some details about the exact solution in Regime III and I on which we will focus. 
In section III, we will concentrate on the computation of the Renyi entropy for regime III, corre- 
sponding to unitary minimal models and in section IV on regime I corresponding to parafermionic 
CFT. To better elucidate the meaning of our formulae, we will conclude the analysis with the 
specific examples of the Ising and 3-state Potts model in section V. Finally in section VI we will 
discuss our results and their meaning. We collect some useful definitions and identities on elliptic 
functions in A. 

II. THE MODEL 
A. Definition 

We consider the restricted solid-on-solid (RSOS) on a square lattice, first introduced in'''. The 
variables on each node are called "heights" and are integer numbers restricted to the interval: 

l<li<{r- 1) (9) 
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A local constraints is imposed to every configuration 

\k-li\ = l (10) 

for each pair of nearest-neighbor The model belongs to the family of interaction round-a-face 
(IRF) models, introduced by Baxter . Each plaquette is given a Boltzmann weight according 
to the configuration of the 4 sites enclosing it: ^2, ^3, ^4)- Here the four sites around the 

plaquette are counted clockwise from the northwest h,l2,h,U- The model can be exactly solved 

h I2 

w{iMM 

h h 

for a proper choice of the weights W where it appears as a consistent restriction of the solid- 
on-solid (SOS) model and hence with the same Yang-Baxter algebra of an eight-vertex model. 
Weights are thus parameterized in terms of elliptic functions and for the details we refer to the 
original work''^. 

At fixed maximum-height r, the phase-space of the model can be characterized by two param- 
eters p, V. The requirement of real and positive Boltzmann weights gives the constraints 

— l<p<l — T] < V < Sr] 

naturally arranged in four, physically distinct, regimes 

I — l<p<0 T] < V < Si] 

II 0<p<l rj < V < 3r] 
III < p < 1 —f] < V < f] 
rV —1 < p < —1] < V < 7] 

The parameter rj is related to p by 

K{p) 

V = 

r 

where K(p) is the complete elliptic integral with elliptic "nome" p. The regimes I, II and III, IV 
are separated by a line of critical points at p = 0. The parameter v can be considered, roughly 
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FIG. 1: The four regimes in the {p, v) plane. 



speaking, as the spatial anisotropy of the interactions in the model and does not enter in the order 
parameters and the critical behavior. So, for fixed regime we will ignore it. The manifolds of exact 
solution will be simply hnes parameterized by p G (—1, 1). By comparison with (6), we have that 
close to criticality, i.e. when |p| ^ 1, p ~ A. 

B. Exact solution 

The exact solution in consists of three parts. 



1 . introduce the corner transfer matrix (CTM) that, once diagonalized, allows reducing the 2D 
configuration sum into a ID sum already at finite size; 

2. perform the thermodynamic limit by transforming the finite size expressions into series in- 
volving gaussian polynomials and then taking the limit as modular functions; 




FIG. 2: The action of the four CTM generates the full partition function. 
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3. sum up the partial traces (for fixed value of the central height) of the CTM obtaining the full 
partition function. 

Here we will briefly review the first two steps, that are functional to our derivation. For the last 
step, we will use a slightly different approach with respect to the traditional one, formulated in 
terms of the dual variables. 

Comer transfer matrix 

The method of the comer transfer matrix (CTM), introduced by Baxter , allows the exact 
solution of lattice integrable models, computing both the partition function and the one-point cor- 
relation function (e.g. the magnetization). As shown in fig. II B, four operators A, B, C, D are 
introduced. A is the partition function of the system restricted to the first quadrant and with fixed 
boundary conditions on the positive x and y axis. Similarly the other operators B, C, D are defined 
in the other quadrants and it follows that 

(pcTM)t:i" = (ABCD);!:,;;^- ^ z = tipctm (id 

The local height probability (LHP) for the height li at the origin can be written as 

Pi = Prob(/i = /) = Z~^Tt (di.^iABCB) (12) 

The CTM formalism becomes particularly useful in integrable lattice models, where it becomes 
possible to fully diagonalize the operator Pctm, hence computing the exact spectrum and, thus, its 
trace. In the RSOS case, the last two sites m + l,m + 2 = N, determine the boundary conditions 
and we will take them as fixed. Once in the eigenbasis, the corresponding diagonal operator p^iag 
can be decomposed as 

Pdiag = RT (13) 

where both R, T are diagonal, but R is a weight that depends only on the height at the origin li, 
while T takes into account the configuration on the whole line / = {/i, . . . ,lrn, ^m+i, They 
can be summarized in the four regimes as follows 
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Regime 


t 


In X 


Ri,i 




II 


2 - r 


47r2 


^(2-t)(2«i-r)Vl6r 




III 


2 


r In IpI 


I 


2 -r 


27r2 


2,l/4+(2-t)(2/i-r)2/8r E{x''\ X~''^'^) 




IV 


2 


r In IpI 



where x) is the elliptic function defined in (A6) and we introduced the CTM Hamiltonians: 

0[/] = f:j^^±^, (14) 

m 

J=l 

(15) 

where 9{x) is the step-function with 9{x < 0) = and 9{x > 0) = 1. 
C. Groundstate structure and critical points 

The two functions (14,15), can be considered as Hamiltonians related to the CTM. In fact, they 
appear as energies for the Id configurations in the trace sum of (11). We can therefore use such 
expressions to deduce the form of the groundstate in each regime, as the configurations having the 
maximum contribution in the trace: since each 2d groundstate is invariant under a southwest to 
northeast translation, it will be enough to fix it on a line li, . . . Moreover, the constraint (10) 
naturally divides the system into two sub-lattices, one with even heights and one with odd ones. 
By specifying boundary conditions (at infinity and at the origin, in a consistent way), we assign 
a given parity to each sub-lattice. But a translation of the whole system by a lattice site gives an 
equivalent configuration, with opposite parity. Thus, we can take the central height li and use its 
parity to classify each ground state out of this trivial Z2 degeneracy. 

In the different regimes we have the following structures^: 

I: There is only one groundstate per each parity of li. (Zj, /j+i) = {n,n + 1) and {k, li+i) = 
[n + 1, n) with n = 

^ To avoid additional spurious degeneracies, in this classification we will assume r to be odd for regimes I and IV. 
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II: There are 2r — 4 groundstates (r — 2 for each parity) of the form of ascending and descending 
sequence from 1 to r — 1: e.g. (Zi, h, ■ ■ ■ , Ir, h-i, Ir-, ■ ■ ■ ■, hr-i-, hr-s) = (1, 2, . . . , r — 1, r — 
2, . . . , 2, 1) and all its translated. 

ni: Also in this case we have 2r — 4 groundstates, where all the odd/even sites have the same 
height: hi = X, ki+i = Y with |X - F| = 1. 

rV: As for regime in we have a groundstate for each couple of available nearest-neighbor values 
except for the regime I groundstate values: thus 2r — 4 — 1 — l = 2r — 6 groundstates. 

It is clear that if there is only one groundstate (for each h parity), then we expect the system to be 
"disordered" and this is true in regime I. Indeed, here the order parameter is independent from the 
boundary conditions, within a given parity of sub-lattices. When there is more than one ground 
state (per parity), the system is in an "ordered" phase. 

The critical points can be understood and identified with an appropriate conformal point^^ : 

• I ^ II critical point: the system passes from a disordered to an ordered phase, where p acts 
like a temperature; the critical point has the conformal structure of parafermion. 

• III o IV critical point: here both phases are ordered and the groundstate degeneracy passes 
from r — 2 to r — 3; the critical point has the conformal structure of (r — 1) -unitary minimal 
model. 

D. Thermodynamic limit 

It is easier to approach the critical points starting from region III and I, thus, from now on, we 
will focus just on these regimes. We are interested in the thermodynamic limit of the replicated 
partition function, introduced in (7). The details of the calculation for a = 1 can be found in 
the original work'", thus here we can concentrate only on the main points and the few modifica- 
tions needed. For convenience, we collect some definitions and the relevant properties of elliptic 
functions in A. 
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1. Regime III 

The finite-size partition function for a replicas is easily obtained from (13), by summing over 
the value of the central height 

Z^= x^]" X„(a, 6, c; x^") , (16) 

l<a<r 

where we singled out the boundary conditions as a = Zi, 6 = /m+i, c = lm+2 and 

The thermodynamic limit can be computed exactly once this expression is rewritten in terms of 
gaussian polynomials , resulting in 

X(a, 6, c; q) ^ lim X„(a, 6, c; q) = {q)^' q^^'' Y fa, ^^^illl^; g") , (17) 
where the q-Pochhammer symbol {q)oo is defined in (A3) and 

The partition function can thus be written as 

= lim y [E(x^xO]°X^(a,6,c;x2°) 

l<a<r 

= x^i^'n-J E (a,^±|^;x^") . (19) 

l<a<r ^ ^ 

2. Regime I 

In this regime, using the table in section II B, the finite-size, a-replicated partition function is 
given by 

l<a<r 

where we introduced 

YUa,b,c;q)^ J2 <f^'^ (21) 

As before the thermodynamic limit is computed taking the limit m — 00 in (21). Unlike the 
regime III, here the sum is not convergent, due to the non-zero energy density of the groundstate. 
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Thus, we factor out the diverging contribution (which amounts to an irrelevant redefinition of the 
partition function normalization), obtaining 



lim g-'^y^la, b, c; q) = (q)^' hM) E{q\ q^^) , (22) 
where the boundary conditions enter only in the function fb,ciQ), defined as 

. s . ^biq) 1 < b <n , 

(g) = { - (23) 

1 n<b<r-2 

,1 2<b<n+l 
kb-i{q) ={ ~ ~ (24) 

Er-biq) n+l<n<r — 1 

where n is the integer part of r /2 and we defined 

q^-\l - q^) 



£biq) = 



1-q 



III. ENTROPY AND PARTITION FUNCTION: REGIME III 

Now that we have introduced the model and the replicated partition functions, it is straightfor- 
ward to proceed with the calculation of the Renyi entropy, using (8). However, before we take on 
the full computation, following , we would like to exploit the known relation between the parti- 
tion function with fixed boundary conditions (both at infinity and at the origin) and the characters 
of primary fields in minimal models . This link will drive the expansion of the entropy around 
the critical point, as we will show in section IIIB. Let us discuss this point in some detail. 

A. Characters of the minimal models 

It was noticed in"*-^''' that the quantity in (17) can be identified with a minimal model character. 
This can be shown by simply rewriting (17), using the sum expansion in (A6), yielding 

X(a, 6, c; q) = gK-^)(-'^-i)+fS-A,,. ^J-i)^^) (25) 

which is the character in the minimal model (r, r — 1) of the primary with conformal dimension 

[dr-a{r-l)f - 1 

= — ^h:;-Y) — ' ^^^^ 
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and central charge 



Here the boundary conditions at infinity are accounted just by the combination 

, b+c-1 

a = , 

2 

Therefore, X{a, b, c; q) is equivalent to a generating function of the Verma module degeneracy for 
the representation fixed by the boundary conditions a, b, c. We stress here that this equality is valid 
only at a formal level: indeed, in the one-dimensional configurational sums, the elliptic nome q 
is a measure of the departure from criticality while in the conformal characters q is the modular 
nome related to the geometry on the torus at criticality. Even though they are both usually denoted 
by q, these are two very different objects. 

Similarly, one should not confuse the formal identification of fields in the Renyi entropy ex- 
pansion, with the operator actually responsible for opening the gap. We recall that regime III can 
be described as (6), that is as a lattice deformation of a minimal model by means of a perturbation 
given by e(x) ~ 0i,3(x), which is known to be both thermal and integrable . As we shall see, this 
operator does not appear among the most relevant ones in the Renyi expansion in regime III. 



B. Fixed central height 

We can now compute the Renyi entropy in a sector where the height at the origin is kept fixed. 
As already pointed out in , we stress that this is done at the level of the corner transfer matrix, 
and so, of the reduced density matrix: if such degree of freedom was fixed at the level of the 
Hamiltonian, it would indeed affect the interaction between the two parts of the system. Rather, 
fixing the height at the origin of the CTM corresponds to selecting a sector out of the whole Hilbert 
space of the model, and taking the groundstate within this projection. Thus, we are measuring the 
entanglement of the lowest energy state within this subspace. In general, these states will be a 
superposition of high energy states and thus the calculation of their Renyi entropy can shed some 
light on their properties. 

In approaching the gapless point, the elliptic nome q = x^", in (19), tends to unity. As this is 
not the best parameterization to extract the leading contributions, we perform the customary dual 
transformation of the elliptic nome, granting us an expansion in the original parameter p, which 
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tends to zero at criticality. Using (A8) and (AlO) for one term in the sum of Eq. (16), we obtain 

= ^ r \ ' (28) 

04 (^,P^)p^V2r(r-l) 

where the index a refers to the fixed value of the height at the origin. As stated in' ', the corre- 
spondence with a conformal character, allows to reinterpret this duality, at the very end grounded 
on the Poisson resummation formula, as the invariance of the torus under the modular group. In 
a CFT, every character can then be expressed as a linear combinations of characters of the dual 
theory'' ' 

- e-/^ = J2 Sffx$f {q = e-) , (29) 

t',s' 

where 

Stf = 2, /IIZ(_l)(*+^)(*'+^') sin ( n^] sin ( n'-^] (30) 



r(r — 1) \ r — 1 J \ r 

is the so-called modular matrix. 

To reproduce this result in our setting, we can expand (28) using (Al) and 



nk Q A T 

1^ i^rJ = E E V = ^ ^-^^^^^'^ = ^+2^' + + / + ' ^^^^ 

71=1 fe=l n=l 

where ^^(n) is the sum of the n-th powers of the divisors of n 



+ (32) 

j<i=l j=l 
j-i=n f=n 

To compare the expansion of (28) and (29), it is convenient to use the parameter truly dual to 
the one used in (25), that is 

g = p§ = ^2^^ (33) 

where z/ = (2 — 2Ai 3)~^ = r/4 is the correlation length critical exponent in Regime III''^. 
Collecting everything we obtain 

In Zi^^ = -^\nq + C',a, + A^adrQ^"^ - S^l^.q^^^ + O (^q^^ ) , (34) 

where 

f nd \ /"Ka\ 

7adr = COS I ^— Y jCOSy — j (35) 
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and 

, / 4 , Trd , -Ka\ 

C„j„ = m , = sm sm — . (36) 

\^2r{r-l) r-1 r J 

is the zeroth-order correction, corresponding to the boundary entropy of'\ 

Using (8), we can obtain the expansion for the Renyi entropy in q, while still keeping the central 

height fixed: 



^g2ar(r-l) _Q,g2r(r-l)j ^ Q ^q7(^^ , (37) 



1 — a 

It is well known^ that one can read off the central charge of the model from the coefficient of 
leading term of the entropy, as in (37). Let us remark, however, that it is a pleasant check to 
notice that the standard conformal result, obtained using the replica trick, remains valid also for 
the minimal models, where the twist operator introduced in the computation does not belong to 
the Kac table of the CFT. 

As suggested in , the sub-leading corrections contain information on the operatorial content 
of the theory and their characters. In fact, from (37) and comparing (36) with (30) we recognize, 
consistently with'\ 

C:,, = ln^ff. (38) 

Indeed, the zero-order term is related to the modular matrix between the primary field chosen 
by the boundary condition and the identity, which is giving the dominant contribution. The first 
correction in (37) is coming, as expected, from the most relevant field. Indeed, we see that: 

Qd,a 
_ '^2,2 

and, coherently, from (26) we recognize that the exponent of the correction is A2,2 = ^ri^^i) - 
The identification with the operators of a Virasoro algebra can continue to higher orders, but one 
should notice that the expansion of the logarithm generates additional terms which do not appear 
in the Kac table, such as the second sub-dominant correction in (37), which is just a 2/^2,2- This 
correction is always dominant over the A3 3 = ^r^^crfy • 

It should be noted here, that the parameter q is microscopical in nature and the entropy is 
usually measured as a function of a thermodynamical parameter, such as the correlation length ^. 
From we know that 

^ = ~\nk'{W)=~l^k'{\q\l) ' ^^^^ 
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where 



n=l 



From these expressions we get the expansion: 



^ 64^2 1536^^^ 2949120^6^^^^^'*' ^^^^ 

which should be substituted order by order in (37). At the leading order, this substitution correctly 
fixes the usual normalization in front of the leading logarithm in terms of the central charge Cr 
and the exponents of the corrections as h = 2A in (5). The rest of the terms, however, which 
strictly vanish in the scaling limit, spoil the possibility of reading and reconstructing the operator 
content of the characters appearing in the entropy in any study at finite lattice spacing. This effect 
is completely analogous to the one discussed in' ^' for the XYZ chain. 



a=l 



C. Full entropy 

To calculate the bipartite Renyi entropy of the model in its true ground state we should sum 
over the central height. Using the dual transformation in the full partition function (16), we have 

z^-Y.e°{y.Vp)zi:\ (42) 

where Z^^ is given by (28). We remark that, while zj^^ has a simple interpretation in terms of a 
character, the coefficients in the sum over the central height in (42) do not. In the previous section, 
since a was kept fixed, the value of the coefficient could be absorbed in the normalization of the 
partition function, but now we cannot ignore these contributions anymore. 

For integer values of a this expression can be handled using the infinite sum representation of 
01 in (Alb), giving^ 

\iaZ„ = Cr-lnpH — Inp + ln — = — > In ( 1 — ) (43) 

48n ' ^ 8 ^ ^2r(r - 1) ^ V ^ 7 ^ ' 



+ In2_^p8"('-i) sm I — \ 2_^sm — - } ^(-1)^^^ sm ((2/c + 1 

7=1 ' a=\ ^ Lfc=0 



r / 



^ For arbitrary values of a, we find the infinite product representation of the 6'-functions to be more convenient, 
ahhough completely equivalent. 
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Take, for instance, a = 1: the sum over a can be computed immediately using the orthogonality 
condition 

r-l 

Eiran arrni r 
sm sm = - dn,m (44) 



r r 2 

a=i 



recovering the partition function of the RSOS model with fixed boundary conditions at infinity^^ 
already in the dual formulation, that once expressed in g = pi gives 



2i = xl-^^A——^, (45) 



UT=i (1 



For general values of a, at our knowledge the sum in (42) can not be computed analytically. 
However, it is possible to obtain its expansion order by order close to the critical point introducing 
the coefficients 

Sa{n, fc) = ^si 

The first few terms give 

' ^ "2 



r— 1 

na nan , 2'Ka 
sm" — sm cos , 

a=l 



Cr l + a 

So, = — Ing + ln 

24 a 



2r Txd 
sm 



r — 1 r — i I — a r 



In 



Sa(l,0) 



1 s„(3,0) / 2 '^d 



^4cos^ - 1^ ga-('-i) + C)^g"'-(^-i)^ , (46) 



1 — a 0) 
Some observations about this expression are in order 

• it remains finite as it should, in the « — )• 1 limit due to the properties of Sa-^i{n, k): e.g. 
Si(3, 0) = 0, 0) = |; note that this is a different mechanism w.r.t. (37), where terms 
with and without a at the exponent appear in pairs and together render the Von Neumann 
limit finite; 

• the leading term remains the same as (37) being dictated by the CFT central charge; 

• since Sa{2n, k) = for all integers n, k, every correction coming from the operator A2,2 
and its descendants disappear and the first sub-leading term is now related to the primary 
field of dimension A3 3. We interpret this cancellation as due to the Z2 symmetry 

/ r - Z 

under which the full partition function (42) is invariant, while the most relevant field, being 
identifiable with the order parameter , is indeed odd. Of course, this implies that all odd- 
operators identically vanish in the expansion of the entropy. In any case, as we already 

17 



pointed out, these corrections are not directly ascribable to the operator opening the gap in 
(6), since in general Ai 3 7^ ^2,2^ ^3,3- This will not be the case in regime I: as we are going 
to show in the next section, in the disordered phase the leading correction seems to be given 
by the same operator opening the gap. 



IV. REGIME I 



We can now turn back to regime I and its bipartite entanglement entropy. As we saw in section 
II C, this regime corresponds to a disordered phase, where local expectation values are independent 
from the boundary conditions. Indeed, the structure of the entanglement entropy is different from 
before. As can be seen from (20, 22), the contribution to the partition functions of boundary 
conditions at infinity factorizes out in the term fb,c{(l)- However, this contribution does not cancel 
out in the entanglement entropy (8) and can bring a finite and interesting contribution. In the 
analysis, we will separate the bulk and boundary contribution and consider them separately: 

So. = + S'i''^ . (47) 

Moreover, looking at (22), we notice that, due to the Z2 symmetry, the fixed central height par- 
titions function for a and r — a are equal and additional relations can be established for certain 
values of r and a for their coefficients in (20). 

In approaching the transition toward regime IV the system undergoes a second order phase 
transition described by the parafermionic conformal field theory-'. Here we summarize the main 
features of these conformal points. 



A. Conformal content of parafermions 

The critical point can be described as the coset 

sl{'2)r-2 
U{1) ' 

with central charge 

(48) 

r 

Beyond the conformal one, these theories enjoy an additional 'Lr-2 x symmetry (which is 
actually enlarged to a Wr-2)- This structure allows to reduce the number of allowed anomalous 
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dimensions to a finite set, even for r > 6, i.e. c > 1. These dimensions are determined by the 
charges (Q, Q) under the two Z^._2 symmetries, since each of them is defined modulus r — 2. 
Following^-, we introduce the two indexes 



I = Q-Q, (49) 
m = Q + Q , 

in terms of which the conformal dimension of the most relevant field in each sector can be param- 
eterized as 

< / < r - 2 , 

1(1 + 2) 

<=V^-i(73^' 0<H</, (50) 

/ — m = mod 2 . 

As a matter of fact, each combination {Af^,Af^,) of dimensions for the holomorphic and 
anti-holomorphic part can correspond to more than one primary field. To resolve this degeneracy, 
one need to look into their representation under the VT-algebra . In particular, within the sector 
neutral under the two Z,._2's, i.e. with {Q, Q) = (0, 0), we have the following allowed dimensions 

Al^) = . (51) 

These fields e^, often called energy or thermal fields, are spin-less, that is A^^^ = A^''^ and the 
identity is e^. We recognize that they are degenerate with the parafermionic operators with (/, m) = 
(2A:,0),i.e. A^,,,^ = ^i^^ 

In regime III we showed the exact mapping existing between the partition function at fixed 
boundary conditions and a conformal character, see (25). A similar relation can be established 
regime I as well, but it is less explicit since the mapping is no more one-to-one: we refer to for 
the precise construction. 

In passing, let us point out that the transition between regime I and II can be described as (6) 
where the gap-opening perturbation is due to the most relevant thermal field ei(x). 



B. Fixed central height 

In order to extract the behavior around criticality, we express each term of the sum in (20) using 
the parameter p, through a duality transformation, as we did in (28) and (42). We recall that in this 
regime — 1 < p < 0. However, following the formulae for the partition function and the entropy 
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are to be understood as depending only on the absolute value of p. Thus, in the following we will 
intend the substitution 

P — >\p\ = -P ■ 
The partition function at fixed height at the origin can be written as 

7W — !^ I 

where, for the moment, we dropped the term fbciq), as discussed. 

The computation of the Renyi entropy is quite similar to what we did in section IIIB. Introduc- 
ing again q following (33), with z/P^ = ^2 — 2A[^^^ = 2(.r-2) ^^^^ regime, we obtain 

Si^^ = + ^) log g + - ^ (g^ - af^) , (53) 



where we defined 



27ra 

7a = 1 + cos , (54) 

r 

/2sin^\ 

= log . (55) 



We recognize that, regardless of the boundary conditions, the exponent of the leading correction 
corresponds to A^q in (50), which is not the smallest one. Since this phase is disordered and 
we do not break it explicitly in the computation of the entanglement entropy, we expect only 
neutral fields under the Zr_2 symmetries to enter in (53). Thus, we find it natural to interpret 
the leading correction in (53) as due to the most relevant thermal operator ei, see (51), which, 
coincidentally, is also the gap-opening operator. This interpretation is further corroborated by the 
observation that, expressing (53) in terms of the correlation length ^ ~ g^^^^, the dimension of the 
leading correction becomes 2A^^\ which seems to be due to a spin-less operator. Moreover, as we 
observed in the introduction of section IV, by changing a we can generate only [r/2] independent 
combinations of primary fields and their characters, which coincides with the number of allowed 
thermal operators in (51). 
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C. Full entropy 



We consider the full partition function, that is obtained with the dual transformation of the full 
sum in (20). One gets 



r-l 



a=l 



(56) 



Also in this case, the series expansion of the 9i function (Al) is useful for integer a = n 



pi 



fl 2^"*"" / 2rj 

c^^lnp H — Inp + In — — — } In ( 1 — p{r-2)n 
8 y/r ^ \ 



12n 



+ ln2_^p"{'-2) (-l)J 

j=0 



r-l 



sm 



vra(2j + 1) 



a=l 



oo 

X:(-l)'=(-p)'^sin((2fc + i; 

.k=0 



aTi\ 
r / 



and again the partition function can be reproduced exactly with (44) 



For general a we can expand the Renyi entropy at desired order 

1 



(57) 



^(bulk) 



'24 



a 



.^logr 1 2g^(l,0) 

logg + ^— + log 

2 1 — a r 



1 ^-(3,0)^^^^ _ ^ 



(58) 



1 — a 0) 

where we see that the leading correction comes from the same a[''^ operator as in (53). 

Thus, we see that, unlike for regime III, in the disordered phase the leading correction is less 
sensitive to the boundary condition at the origin and coincide with the scaling dimension of the 
gap-opening field. 



D. The boundary contribution 

Now we turn to the term related to the boundary heights fbdQ)' from its definition in (23), we 
see that for the set of values of b, c that makes it non-trivial, it gives rise to a peculiar set of terms 
appearing in the expansion of the Renyi entropy 



qibc) 



1 



In 



a 



a{r-2) 



eb (X' 



ln6 I (^'~^)^'^ I o( ^ 
6(lng)2 \lng 



(59) 



We see that the boundary contribution modifies the constant term (boundary term' ) and generates 
sub-leading logarithmic corrections. This result may appear surprising: when local quantities are 
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computed as in (12), such term cancels out, as expected since we are in a disordered phase. But 
in the Renyi entropies, being a highly non-local object, also the boundary appears. Of course, 
in considering the entropy of an actual state, one might need to sum over different boundary 
conditions, possibly with different weights, and the boundary contribution to the entropy might 
change significantly. Thus, it might be pointless to try to provide a CFT interpretation of (59), as 
these boundary conditions might not have any conformal counterpart. However, it should be noted 
that, once expressed in terms of the correlation length, the logarithmic corrections in (59) have the 
same form as those predicted in a CFT with a (bulk) marginal field ' . This is peculiar, since 
ljr-2 parafermions present such a marginal field only for certain given values of r. One can check 
that a massless flow in the theory, essentially due to a free boson, is present, for example, for the 
series 



where the central charge takes the form of a free boson plus a unitary minimal model. Since, the 
logarithmic corrections in the entanglement entropy (59) typically appear for every r > 5, as we 
will see in the next section, these terms must have a different origin. 

V. SOME EXAMPLES 

To better elucidate our results, it is instructive to specialize and consider two particular exam- 
ples: the Ising model (c = 1/2) and the 3-state Potts model (c = 4/5). In fact, within the RSOS, 
we have two possible realizations of these models: one in Regime III (respectively r = 4 and 
r = 6) as unitary minimal models and one in Regime I (r = 4 and r = 5) as parafermions. The 
comparison between the two realizations of the same theory can shed some light on the two phases 
and the nature of the corrections. 

A. Ising model 

The Ising model is arguably the simplest CFT, since it consists of only three operators: 1, a, e. 
Characters of the Ising model are known to arise in the study of the CTM's of eight-vertex models'*"^ 
, which has the same Yang-Baxter algebra as the RSOS. Depending on the choice of parameters, 
the 1-D quantum system corresponding to the eight- vertex model is either an anisotropic XY 
model in a transverse magnetic field, or an XY Z chain in zero field. The entanglement entropies 





(60) 
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of both models have been calculated analytically in the thermodynamic, bipartite limit we are also 
considering here. The first has both an ordered and a disordered phase " , while the XYZ model 
presents only the ordered phase '' ''-^ . In the direct parameter x of section II B, the generalized 
partition function in the ordered phase has been observed to be proportional to the character of the 
spin operator, while in the disordered it is a combination of the identity and energy: 

oo 

oc n (1 + ^'") = 2 x-'^''x. {x') , (61) 

n=l 

oo 

oc n (1 + ^'"-') = [Xi {x') + Xe {x')] . (62) 

n=l 

The Kac table for the Ising minimal and parafermionic model can be summarized as 

Ai = Ai,i = Afo = 1 , 
< A. ^ A2,2 = Ap, = ^ , (63) 
A, = Ai,3 = Ajo = I . 

In regime III, as explained in section III A, tuning the boundary conditions at the origin and 
at infinity, we generate each individual character. The entropy at fixed origin height then reflects 
the operator content of the theory under the duality transformation. Thus, looking at the modular 
matrix (30), we see that if we start with the identity or the energy field, the first correction to the 
entropy (37) comes from the most relevant operator, i.e. the spin operator a. However, if the height 
at the origin is set to a = 2, the coefficients of the a contributions vanishes (as the modular matrix 
has zero element for the (a, a) entry) and the most relevant correction is given by the energy, as 
in'^^'^^ 

As we explained in section III C for the general case, the A2,2 field is odd under Z2 and therefore 
disappears in the full entropy (46), and only the energy and identity appear. Of course, the field 
A3 3, which in general would give most relevant correction in (46), does not appear in the Kac 
table of the Ising model and indeed its coefficients are vanishing. 

In regime I, things are a bit different. By direct inspection of (22), we see that for a = 2 (62) 
is realized and both a = 1 and a = 3 give (61). Thus, as we conjectured at the end of section 
IV B, not every combination of operators and their characters appear. In this case it seems that 
fields with the same parity under Z2 appear together. It is then straightforward to see that after 
the modular transformation the character of the spin operator is never generated and the leading 
correction to the entropy is always given by the energy e, both when fixing the central height at 
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any value as in (53) and by summing over it (58) . Indeed, the identity and the energy are the fields 
in (51). 

Moreover, it is easy to check that the entropy contribution (59) due to the boundary condition at 
infinity is always vanishing since for all the allowed values of b, c in (23, 24) we have fbciq) = 1- 

As a final remark, we notice that for the Ising model, through the identities collected in the 
appendix and some manipulations, the partition functions (19, 20) can be written in a relatively 
explicit way. One simplification arises because, due to the Z2 symmetry, it is sufficient to fix the 
boundary conditions to 6 = 1, c = 2. For regime III we have 



^(iii) 



Xi 



n=l 



(64) 



n=l 



and for regime I 



X12 



X 



Xl (X'") + Xe [X'^) 



n(i- 

n=l 

00 

n(i 



X 



X 



'-)' 



(65) 



n=l 



In regime I, the coefficients of the a = 1 and a = 3 terms are equal and opposite and we see that 
for a = 2m — 1 the partition function is simply proportional to the one found in the disordered 
phase of the XY model (62), consistently with the fact that this regime is also disordered. We 
also notice that for a = I, the coefficients in (64, 65) have the same form as (61, 62). Thus, 
the partition function can be formally written as a bilinear in the characters of the model. This 
reminds us of what observed in' ' and we take it as further indication that the character structure of 
the CTM in integrable models is mostly due to the analytical structure that permeates this beautiful 
construction, and not on some underlying Virasoro algebra. 

B. 3-state Potts model 



The operator content of the minimal model in this case is given by: 

.1^12212713 

^ '40' 15' 8' 5' 40' 3' 5' ¥' 



(66) 
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The most relevant is A2,2 = ^, which is the one appearing in (37), except for a = 3. As usual, 
this field cannot enter in the full entropy obtained summing over the central height, and the leading 
correction comes from the next relevant operator with A3 3 = j^. 

In the parafermionic realization, we only have four allowed conformal anomalies 



(53). We notice that in this case, as it was for the Ising model, only two thermal operators (51) 
exist (k = 0, 1) and the boundary corrections again disappears for any allowed choice of b, c. 

We see that, contrary to the Ising example, here the two different (minimal ordered and 
parafermionic disordered) realizations of the 3-state Potts model have different corrections in the 
entanglement entropy. 

The CTM's spectra of the 3- and 5-state Potts model have been calculated numerically in with 
a DMRG approach and an impressive agreement with the analytical expectations was found, also 
in the presence of integrability breaking terms, sufficiently close to criticality. This indicates that 
our results for the entanglement entropy should also remain valid under the same conditions. 

It is worth to recall here that this parafermionic realization of the 3-state Potts model appears as 
the ferromagnetic phase in the Fibonacci chain . As stressed in , due to the topological symmetry 
present in the quantum realization, all the relevant perturbations are forced to vanish, and the 
critical point is topologically protected. It means that, in the RSOS, the topological symmetry is 
restored only at the gapless points. Therefore it would be interesting to compare our predictions 
with the numerical data for the entanglement entropy coming from the anyonic chain, once a 
perturbation breaking the topological charge is turned on. 

VI. CONCLUSIONS 

Following the suggestion put forward in , we show how to calculated the bipartite Renyi en- 
tropy in the thermodynamic limit of a set of models known as RSOS. The method we employed is 
quite general and powerful and requires just the knowledge of the structure of the Corner Transfer 
Matrix eigenvalues of the system under consideration. In our case, the model being exactly solv- 
able, the CTM spectrum is fully known analytically, thanks to''^ However, generic systems close 
to criticality are expected to organize their CTM eigenvalues according to the CFT reached at crit- 




(67) 



The leading correction in the entropies (53, 58) is coming from A^''^ = A2 



pf _ 2 

'2,0 — 5' 



as we found in 
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icality, as seen, for instance, in '^. If one was able to determine the coefficients in the expansion 
of the CTM in terms of characters of the CFT, the approach we used in this work would apply 
directly. 

Beside it feasibility, this study of the RSOS was motivated by the fact that this model provides 
a lattice realization of all minimal and parafermionic conformal models. It is remarkable that a 
single system can realize such a variety of phase transitions and thus its entropy provides a unique 
case study for the approach to criticality in 1 + 1 dimensions. The different CFT's are realized 
by varying an integer parameter r, while the continuous parameter p (or its dual x) measures the 
departure from criticality. Furthermore, the boundary conditions play an important role in fixing 
the phase under consideration. 

We were thus able to compute the dependence of the Renyi entropy on p and to study its 
behavior. The expansion of the entropy in regime III is given by (37), if we project the Hilbert 
space on a subset specified by fixing the central height in the CTM, and by (46) for the absolute 
ground state. In regime I, we have (53) for the projected case and (58) for the whole case, with the 
addition of the boundary term (59), when present. 

Our results confirm the expectation in (5), according to which, approaching the critical points, 
we have a leading logarithmic term with a universal prefactor (set by the conformal anomaly), a 
non universal constant term, and power- law corrections with non-universal coefficients. We related 
the exponents of the corrections to the conformal dimensions of one of the critical fields. The 
leading correction always has the form of an unusual correction, using the terminology of'", and 
its dimension is that of the most relevant field allowed. By changing the boundary conditions on 
the RSOS, we can select different states for which we calculate the bipartite Renyi entropy, and we 
noticed that certain corrections can be suppressed and thus the leading term can be determined by 
different operators. In particular, we found that symmetry considerations prevent the appearance 
of the most relevant field in the Renyi entropy of the the absolute ground state. In the case of the 
minimal models, where the most relevant field ^2,2 is the order parameter, the leading contribution 
is given by the next most relevant operator, that is A3 3. For parafermionic model, the effect is even 
more dramatic, because the Zat symmetry seems to select only certain fields and the first correction 
generally comes from the most relevant operator neutral under the symmetry, that is A^''^ in (5 1). 

In our opinion, these sort of effects due to the boundary conditions could represent an interest- 
ing possibility for numerical studies in this and other models, where the operator content of the 
theory can be in principle read out, by a proper turning of the boundary conditions. 
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In the parafermionic phase, we also observed the emergence of non power-law corrections, of 
the same logarithmic form (log^)^^" expected in the presence of a primary marginal field in the 
theory. These types of terms were already found in'"* and would be in agreement with a naive 
scaling argument applied to the (log£)~^ terms of'", where the expansion is computed for a finite 
interval of length i. However, we already pointed out that these logarithmic corrections are present 
even when the parafermionic theory does not support a marginal field, and thus we should conclude 
that the origin of these terms is not so simple and might be a lattice effect due to non-conformal 
boundary conditions. 

Another possible interpretation is that in general the relation between the corner transfer ma- 
trices and the Virasoro characters is "accidental", in the sense that is purely due to the analytical 
structure of both quantities. Both are elliptic functions: the latter bi-periodic in real space, while 
the former in parameter space (we remind that the elliptic nome q has a different physical interpre- 
tation in the two cases). When expanded close to the critical point, for consistency the CTM has to 
give the correct central charge of the gapless CFT, and this constraints the structure of the elliptic 
series defining the CTM. Since the same constraint applies to the Virasoro characters, this might 
explain why in general one can write the CTM as a sum of characters and why in the RSOS we did 
not find any connection between the dimension of the operator opening the gap and the dimension 
of the most relevant correction in the Renyi entropy. And it might explain why, playing with the 
boundary condition, one can turn on logarithmic corrections with no counterpart in the CFT. 

Finally, let us remark that the original work''^ on the RSOS spent a considerable effort in devel- 
oping advanced mathematical identities (known has generalized Rogers -Ramanuj an identities) to 
access the partition functions of the model. In our calculations, we overcome some difficulties in- 
volved with summing up Gaussian polynomials, by performing first a duality transformation that, 
in our cases, turned a product of Gaussian polynomials into a sum over exponential one, which are 
easy to handle. We do not know how general and applicable this approach is, but it revealed to be 
quite powerful for us. 
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Appendix A: Elliptic functions and q-series 

In this appendix we recall standard definitions and useful identities for elliptic functions and 
g-series, used in the derivations in the text. For a more detailed treatment and for the derivations 
of the various equalities, we refer the reader to one of the standard textbooks on the topic, e.g."'"'^'. 

First of all, the Jacobi Elhptic 9 functions are defined as 

oo 

ei{z, q) = 2 gi^j^2n + l)z] , (Ala) 

n=0 

oo 

e2{z, q) = 2j2 cos[(2n + l)z] , (Alb) 

n=0 

oo 

e^iz, g) = 1 + 2 ^ g"' cos(2nz) , (Ale) 

n=l 

oo 

e^iz, q) = l + 2 ^(-l)"g"' cos{2nz) . (Aid) 



n=l 



Employing the Jacobi triple product identity 

oo 

{xyW)^ {xy-';x^)^ = E ' (^2) 

n=— oo 

where we introduced the q-Pochhammer symbol 

oo 

(«; 9)00 = n (1 - «^') ' (A3) 

A:=0 

{q)oo = {q; g)oo , (A4) 



28 



one can derive the product representations for the 9 functions: 

oo 

ei{z,q) = 2{q^-q^)^q^sinz^[l-2q^''cos{2z) + q*''] , (A5a) 



n=l 

oo 



92{z, q) = 2 g^)^ cos z\[[l + 2q^^ cos(2^) + g^"] , (A5b) 

n=l 

oo 

e,{z, q) = {q'; q') ^ + 2q'-~' cos(2^) + q'^'] , (A5c) 

n=l 

oo 

e,{z, q) = {q'; q') ^ - 2q'^-' cos{2z) + q'^-'] . (A5d) 

n=l 

In the text we also used the function 

oo 

E{z, x) = {z; x)^ {xz-^; x)^ {x; x)^ = , (A6) 



n=— oo 



where the second equality follows again from (A2). 

The duality transformation for 9 functions can be derived using the Poisson summation formula, 
obtaining the so called Jacobi identities. Once we define q and q such that 

g = e'^^, g = e-^, S>r > , (A7) 

they take the form 

d,{z,q) = -i(ir)5e'^^i(r2;,g) (A8a) 

^2(;2,g) = (-ir)^e^^4(rz,g) (A8b) 

9^{z,q) = {-iT)'^e'^d^{rz,q) (A8c) 

e^{z,q) = (-ir)ie^^2(rz,g) (ASd) 

Finally, it is possible to rexpress the function E{z, x) in (A6) by means of the 6 functions 

E (e^'^ q^) = ig-i/V^^i(z, q) , (A9) 

^(-e'",g') =q-^/^&'92{z,q), (AlO) 

Combining these expression with (A8), it is possible to obtain the expression of the partition 
functions in the dual variables. 
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